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ABSTRACT 

We report the discovery of a bow-shock-producing star in the vicinity of the young 
massive star cluster NGC 3603 using archival data of the Spitzer Space Telescope. 
Follow-up optical spectroscopy of this star with Gemini-South led to its classification 
as 06 V. The orientation of the bow shock and the distance to the star (based on its 
spectral type) suggest that the star was expelled from the cluster, while the young age 
of the cluster (^2 Myr) implies that the ejection was caused by a dynamical few-body 
encounter in the cluster's core. The relative position on the sky of the 06 V star and 
a recently discovered 02If*/WN6 star (located on the opposite side of NGC 3603) 
allows us to propose that both objects were ejected from the cluster via the same 
dynamical event - a three-body encounter between a single (06 V) star and a massive 
binary (now the 02If*/WN6 star). If our proposal is correct, then one can "weigh" 
the 02If*/WN6 star using the conservation of the linear momentum. Given a mass 
of the 06 V star of «30 M Q , we found that at the moment of ejection the mass 
of the 02If*/WN6 star was ~175 M@. Moreover, the observed X-ray luminosity of 
the 02 If*/WN6 star (typical of a single star) suggests that the components of this 
originally binary system have merged (e.g., because of encounter hardening). 

Key words: stars: individual: WR42e - stars: kinematics and dynamics - stars: 
massive - open clusters and associations: individual: NGC 3603 



1 INTRODUCTION 

Young massive star clusters lose their massive stellar content 
at the very beginning of their evolution because of dynami- 
cal few-body interactions (Gvaramadze, Kroupa & Pflamm- 
Altenburg 2010b; Fujii & Portegies Zwart 2011; Banerjee, 
Kroupa & Oh 2012). The ejected stars form the population 
of field OB stars (e.g. Gies 1987), whose space velocities 
range from ~10 kms" 1 (the escape velocity from the poten- 
tial well of the parent cluster) to several hundreds of kms -1 
(e.g. Heber et al. 2008). About 20 per cent of high- velocity 
(>30 kms^ 1 ) field OB stars (the so-called runaway stars; 
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Blaauw 1961) are moving supersonically through the am- 
bient interstellar medium (e.g. Huthoff & Kaper 2002) and 
produce bow shocks, which can be detected in the optical 
(Gull & Sofia 1979), infrared (van Buren & McCray 1988), 
radio (Benaglia et al. 2010), and X-ray (Lopez-Santiago et 
al. 2012) wavebands. 

Detection of bow shocks around massive star clusters al- 
lows us to reveal OB stars running away from these clusters 
(Gvaramadze & Bomans 2008; Gvaramadze et al. 2011c). 
Follow-up spectroscopy of stars selected in this way and the 
geometry of their bow shocks provide a powerful tool for 
linking the newly identified massive stars to their parent 
clusters even in those cases when the proper motion mea- 
surements for the stars are unavailable or unreliable (Gvara- 
madze & Bomans 2008; Gvaramadze et al. 2010a,b, 2011c; 
Gvaramadze, Pfiamm-Altenburg & Kroupa 2011a). This in 
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Figure 1. Left: MIPS 24 fim image of the bow shock detected near NGC3603. The position of the associated 06 V star 
(2MASS J11171292— 6120085) is marked by a circle. Right: MIPS 24 /im image of the highly saturated Hn region associated with 
the star cluster NGC3603 (indicated by a diamond). The positions of the bow-shock-producing star and the 02IP/WN6 star (WR42e) 
are marked by circles (see text for details). At the distance of NGC3603 of 7.6 kpc, 30 arcsec and 5 arcmin correspond to sal.l and 10.9 
pc, respectively. 



turn provides useful constraints on the modelling of dynami- 
cal evolution of young star clusters (e.g. Kroupa 2008; Porte- 
gies Zwart, McMillan & Gieles 2010) and has important con- 
sequences for understanding the origin of the field OB stars 
(Gvaramadze et al. 2010b, 2012a). Bow shocks can also be 
used to constrain mass-loss rates of their associated stars 
(Gull & Sofia 1979; Kobulnicky, Gilbert & Kiminki 2010; 
Gvaramadze, Langer & Mackey 2012b) and could serve as 
probes of the Galactic magnetic field (Gvaramadze et al. 
2011b,c). 

In this Letter, we report the discovery of a bow-shock- 
producing star in the vicinity of the star cluster NGC 3603 
(Section^). Follow-up optical spectroscopy of this star led to 
its classification as 06 V (Section[3]). In Section[4] we show 
that the relative position of the 06 V star and a recently 
discovered 02 If*/WN6 star on the sky is consistent with the 
possibility that both objects were ejected from NGC 3603 via 
the same dynamical event - a three-body encounter. 



2 SEARCH FOR BOW SHOCKS AROUND 
NGC 3603 

2.1 NGC 3603 

NGC 3603 is a very young (~2 Myr; Kudryavtseva et al. 
2012) and one of the most massive (~10 4 Mq; Harayama, 
Eisenhauer & Martins 2008) star clusters known in the Milky 
Way. It is located in the Carina arm at a distance of d=7.6 
kpc (Melena et al. 2008). The cluster contains numerous O- 
type stars and three very massive WN-type stars (Moffat, 
Drissen & Shara 1994; Schnurr et al. 2008), whose initial 
masses could be as large as «140— 170 Mq (Crowther et al. 
2010). All three WN-type stars are located in the cluster's 
core and two of them are short-period («4— 9 day) binary 
systems (Schnurr et al. 2008). 

The existence of very massive (binary) stars in the clus- 
ter and compactness of its core (the radius of the core is 
~0.2 pc; Harayama et al. 2008) suggest that, at least in the 



recent past, NGC 3603 was effective in producing massive 
runaway stars (e.g. Gvaramadze et al. 2009; Gvaramadze & 
Gualandris 2011). This makes attractive a search for OB 
stars running away from this cluster. However, given the 
large distance to NGC 3603, it is unlikely that the existing 
astrometric catalogues would provide reliable proper motion 
measurements for stars ejected from the cluster, unless their 
peculiar velocities exceed several hundreds of kms -1 . The 
runaway status of stars around NGC 3603 can, in principle, 
be revealed through detection of their high peculiar radial 
velocities, but in the absence of reliable proper motion mea- 
surements for these stars it would be impossible to unam- 
biguously prove their relationship to the cluster. This makes 
bow shocks the most important signature of massive stars 
running away from NGC 3603. 



2.2 Bow-shock-producing star 
2MASS J11171292 6120085 

To search for bow shocks around NGC 3603, we used archival 
Spitzer Space Telescope 24 /j,m image of NGC 3603 and its 
surroundings (Program Id.: 41024, PI: L. Townsley) ob- 
tained with the Multiband Imaging Photometer for Spitzer 
(MIPS; Rieke et al. 2004). This image represents an ~ 
0?5 x 2?0 strip extended almost in the eastwest direction 
and centred at «20 arcmin to the west of NGC 3603. Visual 
inspection of the image led to the discovery of five clear arc- 
like structures, which we interpret as bow shocks (cf. Gvara- 
madze & Bomans 2008; Gvaramadze et al. 2011c). However, 
only one of them is opened towards NGC 3603 (see Fig.[T]) 
and therefore could be generated by a massive star expelled 
from this cluster. The stand-off distance of this bow shock 
is ~7 arcsec (or «0.25 pc). 

Subsequent examination of the Digitized Sky Survey II 
(DSS-II; McLean et al. 2000) red band image showed that 
the bow shock is apparently generated by a star with co- 
ordinates: Q J 2ooo=ll h 17 m 12 i :93,5j2ooo=-61 20'08'.'6, which 
is located at «0.262 degree (or ~34 pc in projection) from 
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Figure 2. Normalized spectrum of 2MASS J11171292— 6120085 with principal lines and most prominent DIBs indicated. 



the centre of NGC 3603. This star has a visual magnitude 
of f»15-16 (Zacharias et al. 2004; Lasker et al. 2008) and 
2MASS (Two-Micron All Sky Survey) J, H, and K s mag- 
nitudes of 11.79±0.03, 11.22±0.03, and 10.92±0.02, respec- 
tively (Cutri et al. 2003). In what follows, we will use for 
this star its 2MASS name - 2MASS J11171292-6120085, or 
J1117-6120, in short. 



3 SPECTRAL TYPE OF J1117 6120 

3.1 Spectroscopic observations and data reduction 

To determine the spectral type of Jl 117— 6120 and thereby 
to constrain its distance, we observed this star within the 
framework of our programme of spectroscopic follow-up of 
candidate massive stars revealed via detection of their bow 
shocks (e.g. Gvaramadze et al. 2011c). For observation we 
used the Poor Weather time at Gemini-South under the 
program-ID GS-2011A-Q-88. The spectroscopic follow-up 
was performed with the GMOS-S (Gemini Multi-Object 
Spectrograph South) in a long-slit mode, which provides 
coverage from 3800 to 6750 A with a resolving power of 
~3000. The spectrum of J1117-6120 was collected on 2011 
February 20 under fairly good conditions: with only some 
clouds and a ~1 arcsec seeing. The aimed signal-to- noise ra- 
tio (S/N) of «150 was achieved with a total exposure time 
of 3x900 seconds. 

The bias subtraction, flat-fielding, wavelength calibra- 
tion and sky subtraction were executed with the gmos pack- 
age in the gemini library of the iraiQ software. In order to 
fill the gaps in GMOS-S's CCD, the observation was divided 



1 I RAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory (NOAO), which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. (AURA) under cooperative agreement with the 
National Science Foundation (NSF). 



into three exposures obtained with a different central wave- 
length, i.e. with a 5 A shift between each exposure. The ex- 
tracted spectrum was obtained by averaging the individual 
exposures, using a sigma clipping algorithm to eliminate the 
effects of cosmic rays. The average wavelength resolution is 
^0.46 A pixel" 1 (FWHM r;4.09 A), and the accuracy of the 
wavelength calibration estimated by measuring the wave- 
length of ten lamp emission lines is 0.061 A. A spectrum 
of the white dwarf LTT3218 was used for flux calibration 
and removing the instrument response. Unfortunately, due 
to the weather conditions, any absolute measurement of the 
flux is not possible. 



3.2 Spectral classification of J1117 6120 

Fig. [5] presents the normalized spectrum of J1117— 6120 in 
the AA4000 - 6750 A region, where most of the identified 
lines are marked. The spectrum is dominated by absorp- 
tion lines of Hi and Hell, which is typical of O-type stars. 
He i AA5876,6678, O in A5592 and Cm A5696 absorptions 
are also clearly visible in the spectrum. The relatively high 
interstellar extinction towards the star (« 6 mag; see Sec- 
tion l4.1[) is manifested in numerous diffuse interstellar bands 
(DIBs). A blend of strong Nal A5890, 96 absorption lines is 
of interstellar origin as well. 

The poor observational conditions and the extinction 
towards J1117— 6120 degraded the blue part of the spec- 
trum, which precludes us from using traditional classifica- 
tion criteria. Instead, we utilized the classification scheme for 
dwarf O stars in the yellow-green (AA4800 - 5420 A) pro- 
posed by Kerton, Ballantyne & Martin (1999). Using their 
equation (3) and the observed equivalent width (EW) of the 
Hen A5412 line of 0.99±0.06 A, we found that the star is 
of «06 V type, although the uncertainty of the calibration 
leaves the possibility that it is either of 05 V or 06.5 V type. 
Alternatively, the spectral type of J1117— 6120 can be de- 
rived using the EW(H7)— absolute magnitude calibration by 
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Balona & Crampton (1974), which for EW(H7)=2.28±0.06 
A gives the same spectral type of 06 V. 



4 DISCUSSION 

4.1 Parent cluster of J1117 6120 

To estimate the distance to J1117— 6120, one can use the ob- 
served photometry of this star and the synthetic photometry 
of Galactic O stars by Martins & Plez (2006). The existing 
catalogues (available through the VizieR catalogue access 
tool) provide quite different measurements of optical magni- 
tudes, e.g. 5=17.3 mag (USNO-A2.0; Monet 1998), V=15.7 
mag (GSC2.3; Lasker et al. 2008), and 5=16.2 mag and 
V=15.0 mag (NOMAD; Zacharias et al. 2004). In the lack 
of quality optical photometry, we will use only the 2MASS 
photometry (see Section [2T2)) . whose J and K s magnitudes 
are consistent within the margins of error with those from 
the DENIS (Deep Near Infrared Survey of the Southern Sky) 
database (DENIS Consortium 2005). 

Using the absolute if s -band magnitude and the intrinsic 
J—K s colour typical of 06 V stars of —4.13 mag and —0.21 
mag, respectively, we calculated the _K" s -band extinction of 
Ak _=0.71±0.02 mag and the distance to the star of 7.4!°;;? 
kpcl, which agrees well with the distance to NGC 3603 of 
7.6 kpc. The extinction towards J1117— 6120 can also be 
estimated by matching the dereddened spectral slope of this 
star with those of stars of similar effective temperature. In 
doing so, we found the colour excess E(B — V)=2.07±0.05 
mag. Then using the extinction law from Rieke & Lebofsky 
(1985) and the standard total-to-selective absorption ratio 
Rv=3.1, we found A v =6.42±0.16 mag and ^k 3 =0.72±0.02 
mag. The latter estimate is in excellent agreement with that 
based on the 2MASS photometry. 

Alternatively, the extinction can be estimated by us- 
ing the correlation between the intensity of the DIBs and 
E(B — V) (see Herbig 1995 for a review). Using EWs of 
DIBs at A5780 and A5797 of 1.46±0.04 and 0.44±0.04 A, re- 
spectively, and the relationships given in Herbig (1993), we 
calculated E(B-V) = 2.78±0.32 and 3.08±0.43 mag. These 
estimates somewhat differ from those based on the spectral 
slope and 2MASS photometry. This discrepancy could be 
caused by a foreground region of enhanced number density 
of carriers of the DIBs (cf. Gvaramadze et al. 2012c). 

These estimates along with the orientation of the bow 
shock strongly suggest that NGC 3603 is the parent clus- 
ter of J1117— 6120. Unfortunately, this suggestion cannot be 
proved by proper motion measurements because of the large 
distance to J1117— 6120. Indeed, inspection of the VizieR 
database showed that none of the existing astrometric cat- 
alogues provide significant proper motion measurements for 
this star. 

4.2 J1117 6120: a runaway star from a 
three-body encounter? 

The young age of NGC 3603 implies that J1117-6120 was 
ejected dynamically, either because of a binary-binary or 

2 Note the dominating source of error in the distance is the un- 
certainty in the spectral classification. 



binary-single encounter in the cluster's core. In the first case, 
the most common outcome of the encounter is the exchange 
of the more massive components into a new eccentric binary 
and ejection of the less massive ones with high velocities (in 
general, the trajectories of the ejected stars make an arbi- 
trary angle with each other). In the second case, a single star 
(usually the lowest mass star among the stars participating 
in the encounter) is ejected with a high velocity, while the 
binary system recoils in the opposite direction to the sin- 
gle star. In both cases, the ejected (high-velocity) stars gain 
their kinetic energy at the expense of the increased binding 
energy of the post-encounter binary, which ultimately could 
merge into a single star if its orbit is sufficiently compact. 
Thus, if J1117— 6120 was ejected in the field via a three-body 
encounter then a massive binary or a single merged star 
should exist on the opposite side of NGC 3603 (cf. Gvara- 
madze & Gualandris 2011). 

Interestingly, such a star does indeed exist. This 
02If*/WN6 star, called WR42e, was recently discovered by 
Roman-Lopes (2012). It is located at #i~0.045 degree west- 
northwest of NGC 3603 (see Fig.[T]), just on the opposite 
side of J1117— 6120 (recall that this star is separated from 
NGC 3603 by fl 2 ~0.262 degree). If both stars were ejected 
from NGC 3603 via a three-body encounter, then WR42e 
should be a very massive binary or a single merged star. 
Moreover, one can "weigh" this star using the conservation 
of the linear momentum. Assuming the mass of the 06 V star 
of m«30 Mq (e.g. Martins, Schaerer & Hillier 2005), one has 
that at the time of encounter the mass of the recoiled binary 
(now the 02IP/WN6 star) was M=(0 2 /0i)m«175 M . 

Although the existing proper motion measurements for 
both stars are very unreliabl^E one can constrain their pecu- 
liar velocities using the following arguments. The minimum 
space velocity of a star leaving its parent cluster should ex- 
ceed the escape velocity, v csc , from the cluster's potential 
well, which for NGC 3603 is ^10 kms -1 (here we assumed 
that at the time of the three-body encounter, say 1 Myr 
ago, the radius of the cluster was ~1 pc). Correspondingly, 
the recoil velocity of WR42e is vi^v csc , while the ejection 
velocity of J1117-6120 is U 2 =(02/0i)vi>6O kms -1 . These 
estimates have important consequences for understanding 
the nature of WR42e (see next section). 

The velocity estimate for J1117— 6120 can also be used 
to constrain the number density of the ambient interstellar 
medium, n . Using the wind mass-loss rate and terminal 
velocity typical of 06 V stars (M = 10 -7 Mq yr _1 , Voo = 
2500 kms -1 ; Mokiem et al. 2007), one finds n < 2.6 cm~ 3 . 



4.3 WR42e as a merged binary star 

The ejection velocity of J1117— 6120 should be compared 
with a typical ejection velocity resulting from encounters be- 
tween a very massive binary and a single (less massive) star, 
which for binaries with equal-mass components is «0.8w O rb 
(Hills & Fullerton 1980), where v OI ^ is the orbital veloc- 
ity in the binary. This comparison implies that to expel 



3 Note that the huge proper motion of WR42e given in the PP- 
MXL catalogue (Roscr, Demleitncr Sz Schilbach 2010) is erro- 
neous (S.Roser, personal communication). 
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J1117— 6120 with the velocity of ^60 km s~ , the binary sep- 
aration should be ^15 AU. 

A massive binary of this small separation would be a 
source of strong X-ray emission (e.g. Usov 1992). Following 
Crowther et al. (2010; see their Section 5.1 and references 
therein), one can estimate the expected X-ray luminosity, 
L^ p , of shocked winds in WR42e assuming that this sys- 
tem is composed of two equal-mass components (with the 
mass-loss rates of ~2 x 10 -5 Mq yr _1 and wind velocities 
of 2600 kms -1 ) and that 10 per cent of the shock energy 
contributes to the X-ray emission. We found _L^ tp w5.6 x 
10 34 ergs -1 , which is more than two orders of magnitude 
larger than the observed X-ray luminosity of WR42e of 
L x bs Ri2.3xl0 32 ergs" 1 (Roman-Lopes 2012). Moreover, L$ s 
comprises a fraction of ~2— 4xl0 -8 of the bolometric lumi- 
nosity of the star (see below), which is typical of single stars 
(Chlebowski, Harnden & Sciortino 1989). From this it fol- 
lows that WR42e might be a merger product of the recoiled 
binary system, which is coalesced because of the encounter 
hardening. 

Numerical simulations by Suzuki et al. (2007) showed 
that during the merger process the system loses about 10 per 
cent of its mass, which in the case of WR42e corresponds to 
~20 Mq. The resulting very massive star is a fast rotator 
with higher-than-average helium abundance, which might be 
responsible for the Of*/WN-type spectrum of WR42e (cf. 
Walborn et al. 2010). During the subsequent ~1 Myr the 
star additionally loses about 20—30 Mq in the form of stel- 
lar wind, so that its current mass should be «125— 135 Mq, 
which corresponds to a luminosity of \og(L/ Lq)~6.3— 6.5 
(Crowther et al. 2010; Ekstrom et al. 2012). This value 
should be compared with the luminosity of WR42e. Us- 
ing the 2MASS J and K s magnitudes of this star of 10.18 
and 9.04, respectively, and (J—K s )o=— 0.21 mag (Martins 
& Plez 2006), one finds M Ks =-6.25 mag, which for the K- 
band bolometric correction of — (4.4-^5.2) mag (Crowther 
& Walborn 2011) corresponds to log(L/L Q )^6.2-6.5 (cf. 
Roman-Lopes (2012). 

To conclude, if our proposal on the relationship be- 
tween J1117— 6120 and WR42e is correct then one can ex- 
pect that the peculiar radial velocity of the latter star should 
be about six times smaller than that of the former one. For 
J1117— 6120 we measured a heliocentric radial velocity of 
21.4±6.3 kms -1 , which is an average over the hydrogen and 
Hell lines. After correction for the Galactic differential ro- 
tation and the solar peculiar motion, we found the pecu- 
liar radial velocity of J1117— 6120 of —4.8 kms -1 , i.e. the 
star is moving almost in the plane of sky. [This estimate 
was derived using the Galactic constants -Ro=8.0 kpc and 
00=240 kms -1 (Reid et al. 2009) and the solar peculiar 
motion (Uq,Vq,Wq) = (11.1,12.2,7.3) kms -1 (Schonrich, 
Binney & Dehnen 2010).] Correspondingly, the peculiar ra- 
dial velocity of WR42e should also be almost zero. Radial 
velocity measurements for WR42e are therefore of crucial 
importance for testing our proposal. 
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